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CALCULATION OF EDDY VISCOSITY IN A COMPRESSIBLE TURBULENT 
BOUNDARY LAYER WITH MASS INJECTION AND CHEMICAL REACTION 

Volume I 
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Introduction 


This paper presents the results of recent efforts to calculate 
the eddy viscosity and the effect of mass injection into the boundary 
layer including chemical reactions. 

Although the modified van Driest-Clauser eddy viscosity model 
by Cebeci et al . 1 provides good results for external air flow, the 
adequacy of this model for nozzle flow with . 

~ dllU 

a strongly cooled wall has not yet been verified due to limited 
experimental data. Even for nozzle air flow the modified van Driest- 
Clauser eddy viscosity model does not show the experienced relaminari- 
zaLion tendency in the nozzle convergent section and exhibits an 

unrealistic decrease of eddy viscosity downstream of the throat toward 
the nozzle exit , 2 

In order to include the effects of the past history of flow 
and the strong temperature variation across the boundary layer due to 
chemical reaction or wall cooling, the turbulent kinetic energy 
equation has been introduced for the solution of the Reynolds stress 
in this study. A kinetic energy approach was conducted for incom- 
pressible flow by Rotta 3 , Glushko 1 * , Bradshaw et al. s , Donaldson 6 , 

Malse and McDonald 7 , and Beckwith and Bushnell 8 . Patankar and 
Spalding 9 calculated compressible turbulent boundary layer flow with 
a similar method but could not attain satisfactory agreements of the 

skin friction and the heat transfer results with experimental measure- 
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meats. 


In this analysis, the Reynolds stress and the turbulent kinetic 
energy are related with the gradient of mean velocity according to 
Prandtl and Wieghardt 10 . The turbulent kinetic energy equation, 
derived from the Navier-Stokes equations, is simplified according 
to Rotta 3 * 11 . The constants in the reduced form of the turbulent 
kinetic energy equation are obtained from iterations of numerical 
calculations and comparisons of results with available experimental 
data for incompressible flow with or without mass injection, assuming 
the compressibility effects hidden in these constants to be negligible. 
The profile of the eddy viscosity across boundary layers and boundary 
layer thicknesses calculated for supersonic flow of M = 2.5 are in 
good agreement with experimental data by Squire 12 » 13 . Calculations 
are also performed for subsonic air flow with injection of pure nitrogen 
or a mixture of 4% hydrogen mixed with nitrogen, and for the hydrogen- 
oxygen combustion product flow in a nozzle with hydrogen injection. The 
effects of wall cooling and heating are also investigated for the flow 
without mass injection. 

The solution method was based on the Crank-Nicolson implicit 
finite difference technique 11 * and the chemical reaction was assumed 
to be in local shifting equilibrium. 

Basic Equations 

The compressible turbulent boundary layer equations 2 for steady 

state and two-dimensional and axisymmetric flows are presented in a 
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curvilinear coordinate system, neglecting the transverse-curvature 
effect . The derivation of equations is shown in Appendix A and B. 
Continuity 



Element 

Sa / \ 3a 

p u — M (p v + p’ v* J HI 

7 3y 




(4) 


where j - 0 for two-dimensional flow and j = 1 for axisymmetric flow. 
The eddy viscosity, e , is defined as 

E ~ “ (Pv)' u* / (3u/3yj. ( 5 ) 

The flow is assumed to be calorically perfect and obeys the 
equation of state, 

; _ -P-R T 

M < 6 ) 

where the time mean of the correlation between fluctuating density and 

temperature, p'T’, i s neglected because of its small order of magnitude 
as stated by Harvey et al. 1 ® 
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Turbulent Kinetic Energy Equation 


Considering the continuity and momentum equations for compressible 
turbulent boundary layers, the system of turbulent fluctuation equations 
is written 2 in the cartesian tensor notation as 


3u 


!i 

3x k 


3u, u! 3(p u. 1 il 1 u' 

<pu k ),u i fcT + <pu k ),u i + ITT + — 


3 x, 


9 p ' u' 3 p 1 u' 

— -i + a- 

3*i 


3 x 


j 


+ p* 


3u.' 
3x . 


3 u.' 


9 x . 


u u 


J- L-±-i 

3 X|, 


- 2]x 


3u' 

i 

3x k 3x k 


(7) 


Introducing the previous three equations in Eq. (7) for i= j = 1, 

2, and 3, and defining the turbulent kinetic energy as K = u' u T + 
v’ v’ + wV, the following turbulent kinetic energy equation is 
obtained 2 , using Rotta’s assumptions, 3 , 11 to model the turbulent diffusion 
and dissipation terms. 


_ - 3K 


3K 


3u + 3_ 


pu 3x + (Pv + p'v') 3 y = - 2 (pv)'u' gy 3y 


(lrtapAK I/2 ) f*j 


K 


p K 3 / 2 


- P y 6 — - y 

A 2 A 


( 8 ) 


In Eq, (8) a, 8, and y are constants, and A is the dissipation length. 


In order to relate the Reynolds stress, -(p v) ’ u', with the kinetic 

energy, K, the Prandtl-Wieghardt formula 10 is utilized 
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- (p v) 1 u' 


3u 

3y 


(9) 


= k p A K 1 ^ 2 


The equations can be solved in closed form, when the previous 
three constants, the constant, k, and the dissipation length, A , which 
is a function of the distance from a wall, are known. 

For the ideal (one component) gas, only the equations of continuity, 
momentum, energy, turbulent kinetic energy, and state must be considered, 
and terms including the species gradient can be deleted from Eq. (3). In 
calculating air flow without chemically reactive mass injection the 
Prandtl number is assumed to be constant, and the molecular viscosity is 
considered to follow the Sutherland law: 


y = 2.27 x 10" 8 T i/2 (1 + 198, 6/T) -1 


( 10 ) 


The turbulent Prandtl number is calculated based on the formula by Cebeci : 16 


Pr = o-L_ 

X T 0.44 


[1 - exp (~y t /A + )l 
[1 - exp (-ytp^/ 2 B+) ] 


At the wall, y - 0 S Eq. (11) reduces to 


( 11 ) 



0.4 

0.44 


ll p -1/2 

At r 


( 12 ) 


For a combustible gas mixture the molecular viscosity, y, is cal- 
culated from Wilke’s semi-empirical formula; 17 


» \ a* £ Xj/Xi)- 1 (13) 

i=l j=l 

and the thermal conductivity is obtained according to Mason and Saxena, 18 
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( 14 ) 


* ■ ^2 X ( 1 + 1.065 
i»l 


‘s 

j =1 

J-l 


where the viscosity of each specie, y ± , is taken from Svehla 19 as in 
Reference 14. The laminar Prandtl number is internally obtained from 



(15) 


and the turbulent Prandtl number is calculated as in the ideal gas case 
To simplify the calculation the laminar and turbulent Lewis numbers are 
assumed to be unity. The system of equations (1) through (8) is solved 
including the locally shifting chemical equilibrium calculation 20 to 
relate element and species mass fraction, enthalpy, specific heat, and 


temperature for a given pressure and assigned enthalpies, h ± ° 
The boundary conditions are: 
at the wall, y = 0; 


u (x, 0) = 0 


p v + p' v* = m (x) 

w 

T (x, 0) = T (x) 

w 

ov^y) = • ™ (5 

y=0 w 

K (x, 0) = 0 



W 


fv l 

w e 


w 


(i6) 
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y -*■ 08 ; 


at the outer edge of boundary layer, 

u (x, oo) = U q (x) 

T (x, «) = x (x) 

e 

“« <x> ”> - \ <*> 
K (x, oo) = o 


(17) 


Solution Method 


In reference to experiments by Bradshaw and Ferriss 21 for incompressible 
flow, and assuming the compressibility effect on the dissipation length to be 
negligible, the dissipation length. A, is modeled as 


A/ 6 = 0.205 (y/<5) 3 - 0.586 (y/6) 2 + 0.431 (y/6) (18) 


To avoid occasional negative results of the turbulent kinetic energy K 
in the vicinity of the wall, the boundary layer is subdivided into a region 
very close to the wall and the wake region since the turbulent dissipation 
terms exceed the remaining terms in Eq. (8) due to the combination of 
constants close to the wall. Considering the effect of the molecular 
viscosity in Eq. (8) to be small, the following term can be neglected: 


i- (»**) 
ay \ / 



The eddy viscosity for the wake region results obtained from Eqs. (5) and 
(9) is: 
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(19) 


e 0 - k p A y}/ 2 


3u 



Close to the wall the modified van Driest model 1 based upon the Prandtl 
mixing length theory is applied, and the inner eddy viscosity yields: 


p l 2 


9u 

3y 


( 20 ) 


with the mixing length, 

t = 0.40 y [1-exp (-y/A)] (21) 

The van Driest Damping factor. A," includes the effect of suction or 
mass addition and pressure gradient as indicated by Cebecil 


The inner eddy viscosity, Ej, is used adjacent to the wall until 
the height at which = e. is reached at each axial station. From 
that point to the boundary layer edge* the outer eddy viscosity 
formulation is utilized. 

The three remaining constants in the system of equations are 
assigned the following values: 


a = 0.1/K 

k = 0.6 

Y = 0.36 

After relating a and k by the indicated formulation 2 the constant 
values for k and Y were obtained from numerical calculations, based 
upon a trial and error iteration, which matched 
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experimental data 



by Muzzy 22 for subsonic air flow without mass addition, especially 
the measured eddy viscosity (Fig. 1). Calculation was initiated 150 cm 
upstream of the measurement location to match the given velocity thickness 
90 cm downstream of the computation start point which was considered 
to be the virtual origin of the non- blowing turbulent boundary layer. 

In Muzzy’s experiment 22 the large-scale disturbances were artifically 
created by sandpaper upstream in the test section. 

In order to assess the validity of the present concept, 
calculations were performed for comparisons with Klebanoff’s experi- 
ments. 23 He also used artifical thickening of the turbulent boundary 
layer by covering the first 60 cm of a vertically mounted flat 
plate with sandpaper. Comparison of the measured data and the 
calculated results were made 315 cm downstream from the leading edge. 
According to Klebanoff 23 the turbulent boundary layer virtually 
originates 430 cm upstream of the test station. Thus computation was 

started at this virtual origin. The experimental free stream velocity 
and momentum Reynolds number are U e = 10.67m/sec, and R 0 = 6900, 
respectively. The calculated Reynolds number is R 0 = 6960. The analytically 
calculated eddy viscosities and their corresponding measured values 23 are 
presented in Fig. 2. The agreement with Klebanoff’s data is satisfactory. 

Since Muzzy's experiments 22 include uniform nitrogen gas injection 
from a 60 cm long porous wall, calculations were conducted to obtain 
the relation between the constant, Y , and the non-dimensional ized 
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mass injection rate, F = m w /p e U e . As shown in Figs. 1 and 3, variation of 
the constant y for a given mass injection rate provided good simulations 
of measured values. 

So far subsonic flow was considered. Verifying the assumption 
that compressibility effects on the dissipation length and the con- 
stants are negligible, analytical results were compared with data 12 * 13 
for compressible flow of M = 2.5. As Figs. 4 and 5 indicate, the 
calculated results are in good agreement with measured values by 
Squire. ^ Differences between the calculated and measured eddy 
viscosities non-dimensionalized with the boundary layer thickness 

6 0.995 3re With±h ex P erimenta l tolerances (Fig. 4). In using the 
displacement thickness 6* in the non-dimensionalized eddy viscosity 
term, it becomes evident that the effect of mass injection disappears. 
This, however, is only true for a constant wall temperature, which 
was present in Squire’s case <T W - 295°K) . Different constant wall 
temperatures would indicate a pronounced effect, since the displace- 
ment thickness is strongly affected by the wall temperature. In 
Fig. 5 .calculated boundary layer thicknesses are compared with 
measured values by Squire for a case without mass injection. The 
agreements are remarkably good. The curve shows the profile of the 
turbulent kinetic energy calculated. 
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Effect of Chem ical R eaction in the Subsonic Boundary 
_Layer Due to 4% Hydrogen Injection 


To investigate the effect of combustion on the turbulent intensity, 
the eddy viscosity, and other boundary layer prof iles, hydrogen gas 
diluted by nitrogen gas is uniformly injected from a porous wall into 

the turbulent boundary layer, outside of which dry air is flowing with 

the constant velocitv of . 

' wuuiarige and 

Muzzy, 2 ' and Jonas and Isaacson 25 conductad such experiments, in which 
the hydrogen gas injected is diluted to « by nitrogen gas. The severe 
environment, high temperatures due to strong chemical reaction in the 
boundary layer, prevented them to measure eddy viscosities. In their 
experiments the virtual origin of turbulent boundary layer is far 
upstream of the porous plate and depends on the mass injection rate. 


In order to understand the effects of mass injection and 

combustion on the boundary layer characteristics well and save 

computation time, the present study anaylzes the flow without 

mass injection and the case that the mass,m w = 0.0044g/( S ec. cm 2 ) 

is evenly injected through a 76.2 cm long porous plate, with a 

leading edge. Two modes of mass injection are considered the 

mjectants are 100% nitrogen gas and a mixture of 4% hydrogen and 

96% nitrogen gases by weight. The free stream velocity, U , is 

e 

10*67m/sec s the pressure, - 1 atm, and the temperature, T = 295°K 

e 

Thus, the mass injection ratio if F = m /p U = 0.0035 

wee 



The wall temperature distribution shown in Fig. 6, is almost the same 
as measured by Wooldrige and Muzzy. * For the chemical reaction calculation, 
the following nine species are considered as combustion products: H, H 2 , 

H 2 0, 0, OH, 0 2 , N, NO, and N 2 . 

Calculated results of boundary layer thicknesses are shown in 
Fig. 7. At x 60 cm the 100% nitrogen injection increases the velocity 
thickness, 6, ^approximately 70% of the non-blowing velocity thickness, 
while the combustion due to 4% hydrogen injection causes an additional 
70% increase in velocity thickness compared with the case of 100% 
nitrogen injection. Combustion significantly increases the velocity 
thickness . 

The momentum thicknesses, 0, for the flows with and without 
combustion differ but are in a close range, while the displacement 
thickness, 6*. in the combusting flow is twice as thick as the corresponding 
non-combustion value. The momentum thickness, 0, is reduced due to 
combustion. 

The eddy viscosity profiles are shown in Fig. 8 at three axial 
stations for combusting and non-combusting flows. The magnitudes of 
eddy viscosity for both flows are almost the same. The distributions, 
however, are different. The eddy viscosity profile without mass 
injection at x = 76.2 cm illustrates that both the 
magnitude and distribution are quite different from those of the other 
two cases with mass injection. The peak values of eddy viscosity with 
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mass injection, with or without combustion , is twice as great as 
the value without injection. If the eddy viscosity is shown in a 
non-dimensional coordinate system* such as Figs. 9 and 10, the effect of 
combustion on the eddy viscosity profile is almost concealed. Profiles 
of temperature and velocity are shown in Fig. 11. i n this graph, the 
non-combusting flow has negative gradients of temperature and higher 
velocity gradients at the wall than those of combusting flow. It is 
evident that combustion in the boundary layer plays a similar role as 
injection of mass at the wall, reducing the eddy viscosity close 
to the wall (Figs. 9 and 10). The eddy viscosity close to the wall. 
Fig. 10, is large forn^omb^tlng- caS e, but smHerfor combusting ' 
flow, Fig. 9. This means temperature affects the eddy viscosity. 
Comparing the temperature profiles in Fig. 11 with those measured by 
Wooidrige and Muzzy 2 * and Jones and Isaacson, 25 one can deduce that the 
shifting equilibrium assumption is valid for the combustion of hydrogen 
air. The flame zone shown by their experiments is very thin. 

The difference between the velocity and density profiles in 
a combusting layer and those for a non-combusting layer are illustrated 
in Fig. 12. The second derivative of velocity for the combusting 
flow is positive below and around the flare sheet, while that 
for the non-combusting flow is negative across the boundary layer. 

This result coincides with the experimental velocity profile by Jones 
and Isaacson. 25 Downstream at the axial distance of x = 50.8 cm 
the velocity profiles of the combusting flow tend to collapse 
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towards a single curve, but the density profiles are still non-similar. 

The three velocity curves in Fig. 12 for the cases without injection 
show that heating of the wall is equivalent to injecting mass at the 
wall. 

To understand the great difference of displacement thickness, 

6*, between the combusting and non-combusting flows, mass flow profiles 
are exhibited in Fig. 13. it is obvious that the mass flow profiles 
differ significantly for two cases with injection due mainly to the 
difference of density profiles (Fig. 12). Temperature 
increase due to the combustion is a cause of thickening the displace- 
ment thickness, 6*. The velocity distribution, as shown in Fig. 11, 
is effected by both the distributions of eddy viscosity and density, 
increasing the velocity thickness (Fig. 7). Thus the mass flow 
profile has a dual effect of density, directly and indirectly. 

The curve without mass injection, which is below the profile 

with 100/£ nitrogen injection, shows that the displacement thickness increases 

due to mass injection considering the definition of displacement 

thickness. Effects of mass injection and combustion on boundary 

layer characteristics are summarized at x = 76.2 cm in Table I. 

In Figs. 8, 9, and 10, the magnitudes of eddy viscosity for 
combusting and non-combusting flows are very close while the distributions 
were different. This seems to indicate that combustion does not increase 
the eddy viscosity. The turbulent kinetic energy, however, increases 
significantly due to mass injection and combustion as shown in Fig. 14. 
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This figure indicates that the mass injection increases the turbulent 
intensity by 50% and combustion with 4% hydrogen injection enhances 
it almost three times as much as the peak intensity without injection. 
Thus, it may be concluded that in spite of the increase of turbulent 
kinetic energy due to combustion the decrease of density prohibits the 
increase of eddy viscosity. Considering Eq. (19), the eddy viscosity is 
proportional to the density and the square root of the turbulent kinetic 
energy. Thus, the effect of the density distribution is more 
significant than the turbulent kinetic energy. Comparing curves 

without mass injection in Fig. 14, „ notice that the heated wall decreases 
the turbulent kinetic energy close to the wall. 

The friction profiles at x - 76.2 cm are compared for the 
combusting and non-combusting flows in Fig. 15. The skin friction of 
combusting flow is only IK compared with non-combusting flow, because 
the velocity gradient at the wall of combusting flow is much smaller 
than that of non-combusting flow as shown in Fig. 11. The value of 
skin friction, for flow without injection of course has the largest 
value as shown in Fig. 15. In Table I and Fig. 12, the skin-friction 
V is aiso educed when the wall temperature increased. 

It requires a long distance for the combusting flow to reach an 
almost equilibrium profile. As shown in Fig. 12, velocity profiles were 
almost equilibrium downstream of x - 50.8 cm, but density profiles 
were still non-equilibrium. The distributions of element mass fractions 
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at the wall are shown in Fig. 16. The equilibrium condition at the 

is attained around x — 50 cm. The small change downstream of 

x = 50 cm is due to the wall temperature distribution. 

The profiles of mixture ratio, F/0, are shown in Fig. 17, where 

the mixture ratio is defined as the ratio of hydrogen element to the 

sum of oxygen and nitrogen elements by weight. Fig. 18 shows the 

species mass fraction profiles at x = 76.2 cm. The mass fractions of 

H 2 , H 2 0, and N 2 at the wall are Y r = 0.012, Y H 0 = 0.228, and 

2 2 

= 0.760, respectively. 

Fig. 19 shows the distributions of skin friction coefficient, 

Cj, and skin friction, T w , in flow direction. As known already from 
the velocity profiles, the skin friction is reduced significantly due to 
combustion in the boundary layer. Table I also indicates that the skin 
friction with 100% N 2 injection drops to about one-third of the value 
of isothermal flow without injection. The skin friction of combusting 
flow is less than 10% of the value of the isothermal case without in- 
jection. Even without mass injection, the skin friction becomes smaller 
for higher wall temperatures and increases for cooled walls. The 
effects of cooler walls on velocity profiles are shown in Fig. 12, that is 
the velocity curve indicates that the cold wall has a similar effect on 
the velocity profile as suction at the wall. We have already known that 
a hot wall plays a similar role to mass injection. Referring to Fig. 
again and the velocity thickness shown in Table I, the turbulent kinetic 
energy has a smaller value close to the wall for a hot wall compared 

with a cold wall. The same can be said of eddy viscosities. Fig. 8, 
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where the eddy viscosity without injection and with a cold wall is 
larger close to the wall than for a hot wall. Thus, the velocity gradient 
at the cold wall is steaper than for the hot wall (Fig. 12). As mentioned 

previously the velocity profile is affected through both the eddy viscosity 
and density profile . 

If the weight percentage of hydrogen in the gas mixture injected 
increases to 8% or more using the same mass injection rate as before, 
the boundary layer is blown off as a result of zero velocity gradient 
at the wall at a certain distance from the leading edge. 

Effect or Pu re Hydrogen Injection into 
the Supersonic Boundary Laver 

The flow characteristics in nozzles with divergent half angles 
of 10° and 15° and with uniform hydrogen injection at the wall were 
investigated by Omori. 2 The present study was performed on the boundary 
layer flow in a rocket nozzle with the divergent half angle of 31.5° 
and the variable mass injection rate of pure hydrogen along the 
nozzle contour o 

The thrust chamber geometry, a nozzle built by Pratt & Whitney Aircraft, is 
shown in Fig. 20. The chamber stagnation pressure, P Q , is 2.135 x lO^/m 2 

(3097 Psia); the temperature, T Q = 3680°K, and the mixture ratio of 

oxygen to hydrogen, 0/F = 6.29. The experimental distributions of the 

static pressure and the wall temperature with hydrogen injection 

shown along the nozzle axis in Fig. 21. As shown in Fig. 22, the 
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hydrogen injection rate is uneven in axial direction with the 
highest rate around the nozzle throat. Using these distributions and 
the chamber initial conditions outlined above, calculated results of 
the boundary layer velocity, momentum, and displacement thicknesses 
are illustrated in Figs. 23 and 24 for two cases with and without 
injection, where the same wall temperature distribution in Fig. 21 was 

assumed even in the case without injection. 

For the flow without injection, the velocity boundary layer 

thickness (Fig. 23) decreases in the nozzle convergent section due 

to the reduction of eddy viscosity in flow direction for reasons 

shown in Fig. 25. The velocity boundary layer thickness with hydrogen 

injection in Fig. 23, however, has an irregular distribution because 

of the uneven mass injection. In the nozzle divergent section the 

velocity thickness without injection increases in flow direction 

except for a small disturbance which is caused by the wall temperature 
distribution. It is noted from Fig. 23 that the effects of hydrogen 

injection and upstream history on the velocity boundary layer thickness 

are significant. Fig. 24 indicates that both momentum and displacement 

thicknesses are greatly affected by hydrogen injection. The reason for 

the momentum thickness to become negative with hydrogen injection, 

Fig. 22, is explained below. With injection, the temperature across 
the boundary layer decreases due to ' the” high heat capacity of hydrogen. 

The hydrogen injected reacts with combustion products generating heat. 
This heat energy is consumed to accelerate the flow with an over- 
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shoot in the boundary layer as well as to increase the temperature 
gradient over a certain distance from the wall (Figs. 26, 27, and 28). 

This velocity overshoot shown in Figs. 26 and 28 is caused by the 
combined effects of the favorable pressure gradient of the free-stream 
and the chemical reaction due to hydrogen injection in the boundary 
layer (Jones and Isaacson 25 ). As Fig. 28 shows, the skin friction 
coefficient increases with hydrogen injection and chemical reaction 
in the flow with a favorable pressure gradient. The velocity gradient 
at the wall in Figs. 26 and 28 is greater for the flow with injection 
than without injection. Fig. 27 shows temperature profiles in the 
subsonic region, and Fig. 29 in the supersonic region. 

Profiles of the turbulent kinetic energy are shown in Fig. 30. 

The turbulent kinetic energy without injection has a smaller value 
than that with hydrogen injection except in the immediate vicinity 
of the wall. It is evident from Figs. 29 and 30 that the turbulent 
kinetic energy distribution is related to temperature. Fig. 5 showed 
results of a nozzle air flow with wall temperatures at room condition and lower 
free-stream temperature. The curves in Fig. 30 are the results of the 
accelerated nozzle flow with the wall temperature of 1056°K and the 
free-stream temperature 2376°K. 

Eddy viscosity profiles in the supersonic region are shown in 

Fig. 31. Both the magnitude and the distribution of eddy viscosities 

differ considerably between the flows with hydrogen injection and 

without injection. The magnitude of eddy viscosity increases due 
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hydrogen injection. This result is contrary to a previous report by 

this author, 2 where the mass injection rate was even. The present 

nozzle, however, has a peak injection rate of m = 2. 2g/ (cm 2 * sec) 

w 

around the throat (Fig. 22) and a larger nozzle half angle. This 
uneven injection and the difference of nozzle geometry, which causes 
the different favorable pressure gradient in flow direction, are 
believed to increase the eddy viscosity. The eddy viscosity profile 
with hydrogen injection has a minimum point in the middle of the boundary 
layer, where the temperature increases sharply due to the chemical 
reaction of injected hydrogen with combustion products. 
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Conclusions 


Coupling the turbulent kinetic energy equation with boundary layer 
equations, the characteristics of compressible turbulent boundary layer 
flows are solved for cases with or without mass injection and combustion. 
The mechanism of turbulence has been considered in such a way that the 
variation of turbulent production, diffusion, dissipation terms in the 
turbulent kinetic energy equation were examined by changing the constants 
until the results would simulate Muzzy’s experimental data. After 
obtaining the relations for each term in the turbulent kinetic energy 
equation as mentioned above, this modeling of the turbulent kinetic 
energy equation was verified, comparing the calculated results with sub- 
sonic and supersonic experimental data by Klebanoff and Squire, 
respectively. 

To investigate the effect of combustion in subsonic flow, a 
calculation was performed for air flow with mass injection of 4% hydrogen. 
It was found that combustion significantly induces turbulence. The 
magnitude of eddy viscosity in the combusting flow is almost the same 
as that in non-combusting flow with the same mass injection rate of pure 
nitrogen, while the distribution of eddy viscosity differs. This is 
attributed to the difference in temperature profiles between combusting 
and non combusting tlows. Combustion has a similar effect on the velocity 
profile as if the mass injection rate is increased. That is, the 
combusting flow has a smaller velocity gradient at the wall than the 
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non-combusting flow with the same mass injection rate. As a result, the 
skin friction of combusting flow is much smaller than that of non- 
combusting flow, for a case without pressure gradient in flow direction. 

Calculations were also conducted for combustion product flows in a 
rocket thrust chamber with pure hydrogen Injection from the wall simulating 
transpiration cooling. The injected hydrogen decreases the temperature in 
the boundary layer due to its high heat capacity. When the mass in- 
jection rate is varying in axial direction, 

the eddy viscosity becomes larger than for the case without injection 
although nozzle geometry has an affect also. The distortion of eddy 
viscosity occurs, when the temperature increases sharply due to the 
chemical reaction of injected hydrogen with combustion products in the 
boundary layer. The velocity gradient at the wall for the case with pure 
hydrogen injection and chemical reaction is larger than that without mass 
injection for a flow with favorable pressure gradient as prevalent in a 
nozzle. Thus, the skin friction of the flow with hydrogen injection is 
larger than that without injection, when a favorable pressure gradient 
exists. This result is contrary to constant free stream flow conditions. 
Therefore, the pressure gradient in flow direction exerts a significant 
influence on the boundary layer characteristics. 

The effects of wall temperature on the boundary characteristics were 
also examined for a flow without wall mass injection. When the free 
stream velocity is constant in flow direction and subsonic, heating of 
the wall is equivalent to injecting mass at the wall, and cooling the wall 
has the opposite effect. 
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SYMBOLS 


A van Driest’ s Damping factor 

B Mass transfer number: 2F/Cf 

Cf Skin friction coefficient 

Cp Specific heat at constant pressure 

D Diffusion coefficient 

F Mass injection ratio, to /p e U p 

w 

H Total enthalpy, H - H + H* - h+ u 2 /2 


Enthalpy, h=h+h’ =VhY 

i i 
i=l 


Enthalpy of species i, = h ± + h’ = f C p4 dT + h^ 


Turbulent kinetic energy per unit density, K = u’u' + v'v' + 7 i w r 


Constant 


Le Lewis number 


i. Prandtl mixing length 

M Mean molecular weight 

<1 

m w Mass addition rate from the wall 

n Number of species 


P , p Pressure 


P c Chamber stagnation pressure 


Prandtl number 


S* Heat transfer rate into the wall 

R Universal gas constant 

r mi Msss ratio of element m in species i 


tty Nozzle radius 
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St Stanton number 

T Temperature 

u velocity in x direction 
v Velocity in y direction 


Mole fraction of species i 


x 


Y. 

1 


a 


Distance along the contour 
Mass fraction of species i 

Normal distance from the wall 

Axial distance 

Constant 


Oju Element mass fraction, a m = ^ r m ^ 


i=l 


B Constant 


Y Constant 

*5. 990 Velocity thickness, $ = y at u/u £ = 0.990 

^.995 Velocity thickness, 6 = y at u/u e = 0.995 

CO 

6* Displacement thickness, <5* “ / (1 - pu/p e U e )dy 

£ Eddy viscosity 


CO 

(P u /P e u e ) (1 - u/u e )dy 

A Dissipation length U 

^ Thermal conductivity 

V Molecular viscosity 

p Density 

T Shear stress 


Momentum thickness, 0 


-/ 
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Subscripts 


e Boundary layer edge 

I Transpiration coolant 

i Species 

m Element 

T Turbulent 

w Wall 



( ) 1 Fluctuating term 
( ) Time averaged quantity 


Notations in Eq. (11) through Eq, (14) should be consulted with Ref. 14 
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Figure 2. Eddy Viscosity Prof iles^for Klebanoff’s Experiments 
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Figure 4. Eddy Viscosity Profiles for Squire's Experiments 

- 33 - 












Figure 



Air Flow with Combustion due to Injection of a Mixture of Hydrogen and Nitrogen 
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Figure 10. Non-Diraensionalized Eddy Viscosity Profiles without Combustion 
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Figure 13. Mass Flow Profiles 
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Figure 14. Turbulent Kinetic Energy Profiles 
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Figure 17. Mixture Ratio Profiles 
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Figure 23. Calculated Distributions of Velocity Thickness 
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Figure 24. Calculated Distributions of Momentum and Displacement Thicknesses 






Figure 25 « Comparison of Eddy Viscosities with and without Hydrogen Injection 
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Figure 26. Velocity Profiles with Hydrogen Injection 
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Figure 29. Comparison of Temperature Profiles with and without Hydrogen Injection at H 
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Figure 30. Comparison of Turbulent Kinetic Energy Profiles with and without Hydrogen at M=2.35 and 2.8 
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Figure 31. Comparison of Eddy Viscosity Profiles with and without 


Hydrogen Injection at M=2.35 and 2.8 



Table I EFFECTS OF MASS INJECTION AND COMBUSTION 

ON BOUNDARY LAYER CHARACTERISTICS AT x = 76.2 cm 



COLD WALL : WALL HAS A CONSTANT TEMPERATURE OF 150 °K 

ISOTHERMAL: WALL TEMPERATURE IS SAME AS MAIN FLOW TEMPERATURE 

HOT WALL: WALL TEMPERATURE DISTRIBUTION IS AS SHOWN IN Fig. 6 









APPENDIX A 


TURBULENT BOUNDARY LAYER EQUATIONS 


1* Continuity Equation 

The overall continuity equation for steady flow reads 
3 (p u i ) 

— r - 0 (A-i) 

9 xj 

in cartesian tensor notation. 

In turbulent flows this equation must be true at any instant, 
but also on the average. The instantaneous quantities are replaced by 
the time average plus the fluctuating as follows: 


P = p + p' 

an< * u i = u i + ui’ 

Substitute these expressions in Eq. (A-l) , then 

3 t (P + P') (Z ± + Ui ’) ] 

9 

Since 


(A-2) 
(A- 3) 


(A-4) 


(p + P 1 ) (ui -J- Ui 1 ) - p + p Ui ' + p ? U;L + p* Uj/ (A-5) 
taking the time average of this equation and applying Reynolds' rule* 
of averages, we obtain 

(P + P') (u ± + up = P u ± + p Ui + p* + p' Ui ' 

P u i + P u i' + P* + p’ u^ (A-6) 

- P Uj. -I- p" r “ui 1 ' 

Therefore the continuity equation is written as 

Hughes, W. F., and Brighton, J. A., "Schaum's Outline of Theory and 
Problems of Fluid Dynamics," Schsum Publishing Co., 1964, p. 179. 



9 (p u. + p 1 u. ') 

»*i ■ = 0 

9 


(A-7 ) 


This equation is expressed for two-dimensional or axisymmetric flow in 
a curvilinear coordinate system as follows, in which s is the wetted 
length along the wall and y is measured normal to it* 


9 . 5 

3 S I (P u + p' u') [ (p v + p T v') r^] = 0 


(A-8) 


where 


j = 0 for two-dimensional planar flow 
and j = 1 for axisymmetric flow. 

Assuming 


9 (p r u’ ri) 3 (p' v V r J) 

« w 


9 s 


9y 


(A-9) 


Eq. (A-8) becomes 


9 (p u rJ) 3 [(p v + p’ v’) r^ 


= 0 


(A— 10) 


3s 3y 

The velocities, u and v, are those components along s and y, respectively. 

2. Momentum Equation 

The momentum equations are written in a cartesian notation as 
follows** for steady flow: 


3u-i 9 P-? * 

°‘i^ ' -3? 


where the stress tensor, , is 


(A-ll) 


li3^ < A ' 12 > 


*Schlichting, H, , ’’Boundary-Layer Theory", McGraw-Hill Book Co., 1968, p. 223 
**Landau, L.D., and Lifshitz, E. M. , "Fluid Mechanics," Pergamon Press 
1959, p. 48. 
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The quantities y and £ are called the molecular (or dynamic) viscosity 
and the bulk viscosity, and are functions of pressure and temperature. 
Neglecting the bulk viscosity tern in Eq. (A-ll) according to the Shvab- 

Zeldovich f emulation 1 *, we can write Eq. (A-ll) in cartesian coordinates 
as 


Pu 


1 



+ pu 2 




and 



(A-13) 


9u„ 9u 0 

p u — L + p u, — - 
^ 9x^ 9 x 2 



(A-14) 




2 / 8u l 9 u 2 \ 

3 \9xj_ + 9 x 2 / 


Because of the boundary layer assumption +t Eqs. (A-13) and (A-14) reduce 


to 


P U, 



+ 


9u^ 

p u 2 5^ 


9p 9 

9x^ 9 x 2 



(A-15) 


and 


aT =0 (A-i6) 

2 

multiplying U;L by Eq. (A-l) , we obtain 

''Williams, F. A. » "Combustion Theory," Addison-Wesley Publishing Co., 
Inc., 1965, p. 9. 

^Schlichting, H., "Boundary-Layer Theory," McGraw-Hill Book Co., 1968, 
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3 (p Uj) 


(A-17) 


u. 


3 (p u 2 ) 

— — + u- — <* 0 

9x i 9 x 2 


Adding Eqs. (A-15) and (A— 17) gives 

9 ( p u l 3 (p ui u 2 ) 

3 x-^ 3 x 2 

Substituting the following relations into the above equation, 

P = p + p’ 


_ 3P + 9 

3x^ TZq l 3xj 


(A- 18) 


u i " u i + u i* (A-19) 

P = P + P r 

and taking the time averages, we obtain 

3 {(p + p') (u 1 .+. Ul t ) 3 | (p+ p*) (G 2 + u^)} 

3 x, a „ 


3 (p + p ' ) 3 

3 i + 3 x 2 


3 (u x + u^) 


3 X' 


(A- 20) 


According to the Reynolds rule of averages, 


(p + p') (uj + u^') (u-j_ + u^’) 


(p + p') (u^ + 2 u^ U2 + u^ u'^) 


* P u 1 ^ + 2 p u 1 + p u^ 

(A-21) 


+ p' u-j^ u 1 + 2 p’ + p' 


- pu lU;L + pu^ + 2 p’u^ + p’u^u’ 
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(P + P') (U^ u ! 1 ) (u 2 + u’ 2 ) - <p + p'Hu^ + u^ + u’^2 + u'^) 


p u x u 2 + pu lU ' 2 + pu\u 2 + pu'^ 


+ p’ u x S 2 + p’ ^ .+ p'u' 1 u 2 + p'Uj* u£ 


" P«i«2 + P u i u 2 + P' u 2 5 i + P ,u i« 2 + p’u^ ul. 


P + P’ = p 


U. + u. 8 = u_ 


111 

Thus Eq. (A-20) reduces to 


9 (p u u + p u 8 u ? + 2p' u ’ u + p 8 u 7 u M 
1 1 1 1 11 1 l 


9 (P u ! u 2 + pu l u 2 + P#u 2 U 1 + P" u , ? u + P* u 8 u * 
— 12 12 


(A-25) 


Substituting Eq. (A-7) into the above equation, then 


(P Uj. + p' Ul ’) 1 + (p u 2 + p f u 8 ) ^1 


9 / 9u]\ 9 (p u x T u x 8 + p’u^ Ul ' + p'uj 8 uj ) 


9 p 9 / 9 u^ 

9x^ 9 x 2 \ 9x 2 
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(A-26) 


3 (p u 2 ' + p* u 2 + p’u^ u 2 ) 


Since 


and 


p Ul ' u{ + p» u L ' u x + p’uj/ Ul ’ « (p Ul ) f Ul r 


p u * u 1 


1 ~2 ■' P U l’ U 2 + P 'V V = (pU ">' u ’’ 


2 1 


Eq. (A-26) is written in a cleaner form as 

9u 


— - ou- 3ui 

(p u 1 .+ p' Uj/) _1 + (p u 2 + p’u 9 ’) ^ 

1 


9xi ' 2 r “2 3 x 2 


(A-27) 


(A-28) 


(A-29) 


3 p 

3 

f 3U 1 


> 3 

{ ) 

3xj 

+ 9^ 

r ^2 

- (pu 2 ) ' 

r 9x i 

|(pui)’ Ul ' j 


This equation is expressed in a curvilinear coordinate system for two- 


dimensional or ax i symmetric flow,* when the transverse*“curvature effect 


Assuming 


and 


as 


— t \ 

u >3l + 

(p v - 

iP . ± \ 

i 3u 

3s 3y j 

3y 

, , 3u 

p u 3l 

< < 


3u 

9y 


" (pv) ’ u ' ‘al 


3u 


j (pu) ' u’ | (A-30) 


3y 


(A-31) 


~ir { (pu)l u '} < < { <pv> 


(A-32) 


*Schlichting, H. , "Boundary-Layer Theory," McGraw-Hill Book Co., 1968, 
p. 223. 


i 


- 67 - 



I 

I 

we finally obtain the momentum equation in s— direction. 


3u. 

p u 3s + v + 





9u _ 

3y ~ (pV)’u') 


(A-33) 


The y- direction momentum equation (A-16) yields 


p + p’ = constant 

across the boundary layer at each station. 
Energy Equation 


(A-34) 


3 * 


The general energy equation *,** for steady state is written 
for two-dimensional or axisymmetric flow in cartesian tensor coordinates 


as 


P ui 


3H 


+ p u 2 



9 Hj 3H 

9x 2[ p r 9x 2 


+]iCl - — ) u 

P r 9 x 2 


where 


+ P_ “1) 

* r 


S hi 

i=l 



H 

h 



i=i Vi 


(A-35) 


hi 


T 

/ 

To 


C 

Pi 


dT + h.° 


x 


(A-36) 


Hughes, W.F., and Brighton, J.A., "Schaum's Outline of Theory and 
Problems of Fluid Dynamics," Schaura Publishing Co., p 196 

Lees, L., "Convective Heat Transfer with Mass Addition and 
Chemical Reactions," 3rd AGARD, 1950, p457. 
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p r P "^pA 

?P = Jl Cp i Yi 

L e = PDTIp/ ^ 

u * y (T) 

A = X (T) 

Adding Eq. (A-l) multiplied by the total enthalpy, H, and Eq. 
then we obtain 


9(p Ul H) 9 (Pu 2 H) 

9xi 9x2 


3 

9x2 


MS- +U (i-^) 


9x2 


P. 


9u^ 
*1 3x2 


+ 


fc a '- 1> 


n 

2 

i=l 


hi 



Substitute the following relations 

Uj = Uj + u j ! 

P P + P* 

H i + H' 

Y = Yi + Y i’ 

hi = ^ + hi’ 

into Eq. (A-37), and take the time average, then 


(A-36) 
(A-35) , 


(A-37) 
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3[(p + p'Huj + uj'KH + H* ) j 3ICp- h p’Ku2’+ U2') (H + H’) 


3 x^ 


3 X2 


3x2 


H 3(H + H') _ 

P r 3 x 2 + P (1 - -i — ) (ui + ui f ) 


3 (ui + ui’) 
3 X2 


u n - 3 (Tj+ Y ± ') 

+ f- (L -1) z (h, + V) — J — — 

r e i=l 1 1 3x 0 


(A-39) 


The left-hand side of the above equation is simplified using the 
turbulent continuity equation (A-7) as 


9[(P + P')(ui + ui ? )(H + H 1 ) 3[(p + p’)(u2 + U 2 *)(H + H’)] 

~ + “ " 
^ X 1 9 x2 

9(P u x H + P 1 U]/ H + P ujJ H* + P* H* ui + P’ U3_’ H') 

3 x^ 

3(p u„ I + p* u ’ H + p u ? H’ + fT r H’ ^ ^7?) 
+ ± ± ± 2 2 

3 X 2 

= (P u x + P 5 u*’) |5- + (p^ + p' u 2 ') — 


+ 3(p uj* H* 4- p 1 H* uj + p' ui* H * ) 3(p u 2 * H f + p'H ? u 2 +p , u 2 , H') 

3 *1 + i""x^ 

(A-40) 
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where 


p Ui’H' + p'H* Ul + p* uj.’.H 1 = (pu^'H’ 

p u 2 ’ H ' + pMT‘ U 2 + p» u 2 »h' - (pu 2 ) ? H* 


(A- 41 ) 


Assuming 


and 


3[(PU!)’ H'] 3[(pu 2 )’ H*] 

c< 

a x x a x 2 


P ? « P U-j 


(A- 42 ) 


(A- 43 ) 


Eq. (A-39> is simplified to 
- 9H 


p l U 1 9x x + (p u 2 + P* u 2^ 3x 


3H 


3 x 2 


t !fr (<pu >'■■>♦ t*a 


iE- + 1 _ HulnH \ 

Pr ( 3 x 2 2 a * 2 ( 


+ i^ ( L e-D 


M_ 3 Y ± 3 Y ± * n 

^ + <v ^ >}J 


Since 


H = h + 


U 1 


(A-44) 


* h + h* + 
- . 5 1 2 


tuj + ul') 


& + ~j~ + h’ + u i u i* + 


«1 


H + H’ 


(A-45) 
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Thus 



Neglecting the terms which include three or four correlations, we 
simplify Eq. (A-48) as 

~(pu 2 ) 'H' = - (p'h 1 u 2 + p'u]/ up 2 + P t^'h' + p u x uj/ u 2 ') 

(A-49) 
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The fluctuating term of static enthalpy, h 1 , is expressed as 


— n n 

h' - h - h - Z hiYj. - Z hiYi 
i=l i=l 


n 


n 


1 ^ 1 ( h i +h i')«i +Y i'> - 


n _ _ 

^(hiYi + h±Yi' + hi'Yi + hi'Yi’) 


- Z (h^ + h^'Y, 1 ) 
i=i 


= + h ± , Y 1 + hi'Yi » - hi'Yi') 


(A-50) 


Substitute the above equation into Eq. (A-49) , and neglect third 
correlation terms, then 


-<pu 2 )'H' ZJ u 2 (p'Yi’ ^ + p 'hi'Yi) + p(hi u 2 ’Yi' + u 2 'hi' | ± )) 


(p'ui' Ui U2 + P u i u i' U 2 ') 


Z { (p'hi' u 2 + p .^'hi') Y ± + (p , Y i ' u 2 + p Uj’Yi')^} 


- Uj ( p- Ul ' U 2 + P u i’u2') 


(A-51) 


- 73 - 



where 


p f hi' u 2 + P u 2 ?h i f = - p’^’h^’ 


p’Yi* u 2 + P W = (Pu 2 ) , Y 1 # - P , u 2 , Y 1 ' 


p'V u 2 + P u i' u 2' = (pU 2 )’uj_’ - P ,u i' u 2' 

Eq. (A-51) is rewritten as 


-( P fe 2 )’H’ £ {[(pu2) , h 1 ’- p’u^' ] Y ± + [(pu 2 )’Y i ' 


- ] h ± } - u-l [ (pu 2 )*u 1 l - p , u 1 'u 2 ' ] 


Let us define the turbulent thermal conductivity, 1^., and 
dlffusivity, D,j , as 

B _ g “ 

- [(pu 2 )'h i ' - P , u 2 , h i ' ] Y t = ,A T 


- [ (pu2) 1 Yi' 


p'u2 , Yl' ] 




and assume 


(pu 2 )'u2 ! » p r u^'u 2 ' 


and 


.9_. 9 u i i &2 W 

3x 2 Ul 3 x 2 ^ 2 3x 2 



3Yj \ 9 

3x 2 ' >:> 3x 2 


<V 



(A-52) 

(A-53) 

(A-54) 

(A-55) 

(A-56) 

(A-57) 

XA-58) 

(A-59) 

(A-60) 
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then, Eq. (A-44) becomes 


(P ui + p» + (p u2+ p* u 2 ’) -2S 

5xi 3x2 


[■ x . - 2 s + 


3T - 


,7- 3y1 


3x 2 Pj> 3x 


Xt T“~ +P I>x £ (hi 3 — ) - u (PU2)*U1* 


3x2 


1®1 


X2 


.-J-) «|a + 


+ y (1- ^ ) ui 


*x 2 


(L e -1) £ (h±|Ii)] 

1=1 3x-> 


(A-61) 


In order to replace the temperature, T, in the above equation to the 
static enthalpy, h and h^, and the species mass fraction Y., consider 
the definition of the averaged specific heat, C p . 

5„ - ? Cp Yi = S (|ii) Yi 

i=l 1=1 dT p 


" r 3 (hiYi) _ 3Yi 

" l-l‘ “IT" ■ hi= 3T> 


3h 


n 3Yi 


3T i=i hi 5T 


(A-62) 


Thus 


^ _3T = 3h_ n _&i 

C p 3x 3 x 2 j _J 1 \ 9x 2 


(A-63) 


Taking the time mean of the above equation considering Eq. (A-60) , 
then 


__3T = 1_ 3h_ n - 3Yi 

3x2 Cp '3x2 ~ i=l hi 3x 2 ' 


(A-64) 
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Now the eddy viscosity, e , the turbulent Lewis and Prandtl numbers, 
Le T and Pr T , are defined as follows: 


t ~ ~(pu 2 ) , u 1 ' 



(A-65) 


and 



(A-66) 


(A- 67) 


Considering Eqs. (A-64)through (A-67), Eq. (A-61) is written in a 
curvilinear coordinate system as 



+ (pv + pV) 

3y 


!-(<U- + -£-) IS 

3y I P r P rT 3y 


[Vl(l- - ) + 
Pr 


e(l-i-)] U |E + 
Pr T 3y 


[^( V l) 


p r T 


e T 


•1)] S 

i-1 3y 


! 


where the assumption below was made. 


p'u' 


3H 

3s 


« 


p’v’ 


3H 

3y 


(A-68) 


(A-69) 


4. Element Equation 

Since there is no generation nor disappearance of atoms 
in the system considered, the continuity conservation equation of 
each atom, m, is written as 

div (pn&im) “ 0 (A-70) 


76 - 



where 


Pm = P°<m 


(A-71) 


turn = 'Um.dif + "W (A-72) 

( Dm : Velocity of element m 

w dif : Diffusion velocity of element m 

iu s Mean velocity 

Thus 

div [potm Gum.dif +ai)] = 0 (A-73) 

This equation is rewritten as 

(pu grad) cijn + div (po^ tu^if) div(pu) « 0 (A-74) 

Considering the overall continuity equation, Eq. (A-l), and 
Pick’s Law 

PlUm.dif = - (PD grad) (% (A-75) 


Eq. (A-74) becomes 

(pu grad)^ - div [ (pD grad)a m ] (A-76) 

For two-dimensional or axisymmetric flow, this equation is expressed 
in a curvilinear coordinate with the boundary layer assumption. 

p u + p v (pD-~a) (A-77) 

ds dy dy dy 
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Multiply djj to Eq. (A-l) , then 

3 (pu) 3(pv) 

3s ± “m ddy = 0 

Add the above two equations, and we obtain 


SCpuOjn) (3pva ) 3__ (pD 3a ) 


3s 


m 

3y 


3y 


. m' 


3y 


Substituting Eqs. (A-2) and (A-3), and 


U 


a m -f 


into Eq. (A-79), Eq. (A-79) becomes 


a[ (P+P 8 ) (u-j-u 8 ) (a a! (^{-p') (7-w*) ^-^’)] 

3S 

= iL_ [ (p4-p * )J) SCanrHam 1 ) 

9 y 3y 

Taking the time avetage of this equation yields 

+ (p Trt'PV) T~~ 

3y 


3ain) _ jUpvTSII 

3y 3y 

where the following left hand side terms have been omitted in 

deriving Eq, (A 82) because of their smallnesses compared with 
remaining terms* 


3otm 

P u t — 
q ® 

= |_ (pD 

3y 


(A-79) 


(A-80) 


(A-81) 


(A^8fi) 


the 
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p'u' 


9a 

9s 


« 


pV 


3 cog 

3y 


CA-83) 


and 


3[ (pu) *<%’] 

3s ” 


« 


ai(pv)'V] 

9y 


(A-84) 
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APPENDIX B 


TURBULENT KINETIC ENERGY EQUATION 


The momentum equation for compressible fluid flow reads 


— a(u-rPui’) _ _ 

(P+P ) ^ + (p+p s ) (u k +u k r ) 


9 (u jhij') 
3x k 


_ 3 (P+P 1 ) a (Tik* TikM 


dXi 


3x k 


(B-l) 


where 


tau Tik = Tik + Tik ? 


= y /lHii + iuk_2 ^ 9u£ 

?\ 9«k 3x± 3 ax/ 3xi, 

molecular viscosity bulk viscosity 


(B-2) 


The continuity equation is 


+ KjotolKuk+nk*) _ 

a T 3x k " 0 


thus 

9 (P+P ' ) . afPMfcf Puk' + P’uk + P’uk') 

at + 5 7" = 0 (b-4) 

k 

Subtract from Eq. (B-4) its time averaged equation, that is, 


3p 

at 


+ 


_a(ftuk+ p’u^) 
a^k 


(B-5) 
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and we obtain 


. 9(P-UkV ± PV-i^ P , uk > ~ P'uk') 
St 3xk 


(B-6) 


Multiplying ui' by the above equation (B-6), and adding this and 
Eq. (B-l), we obtain 


— 3ui — 3ui' 3ui 

p 17 + p~ + p F + 


.3 ui* 
p ~ 


+ °i^ 


+ (P Uk+P Uk' + p 1 Uk + p f Uk') ^ + U 1 ) 

3x k 


+ Ui 1 9(p u k’ + P 1 Uk + P r Uk’ ~ P* uk') 

3 *k (B-7) 

= - *jll p,) + 5 ^ik+^ik) 

3x4 3x 

1 k 

Subtract from this Eq. (B-7) the time average equation, then 


0 It * p 'fr + + ( P u, '+p'ut*p'ui l ' s 


3ui 

3xk 


+ (p uk + P uk’ + p r uk + p' uk’)^ _ (p uk’ +p' uk + p' uk') 


*k 


3xk 


+ . 3(P u k' + P' uk + P' u k ' - p' uk') . 3 (P u k 1 + P' uk + P'uk' ~ P’uk') 


u i 


0x k 

3p* 3 T ik' 

3 X . 


- u ± ' 


x * 9 


*k 


(B-8) 
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Write the same equation for the velocity component * , Multiply 

Eq ' (B_8) by V and the e< * uatioa for Uj * by * ± \ and add these two 
equations obtained t then 




+ u j'(P Uk* + P' Uk + p‘ Uk' 1 - p f Uk') — 

9x k 


+ u^(p Ut 1 + p* Ui. + o' 


9u . 


lit ' — r\ f .. f \ J 


K u k > 3x£ 


+ (p Uk+ puk' +p'uk+ p'uk") - - + 2ui'uj ? _^ puk u k~ t ~P u k -p’uk 1 ) 


9xk 


3 U 4 V 


E(pu k ’+p*u k + p’uj^.’) . Ui i 

9 x k 


L , d (P u k ,+ p 1 Uk-S-P v -p 'uk^J 

fl 

3 (puk *+p uk+p ’ uk 5 -p * u k * ) 

L 3% J Ui 1 

Uj — 

J 3*k 




+ Ui ' + uj’ + ui 8 — ik! 

dK l J 3\ 1 3x, 


(B-9) 


After time averaging Eq. (B-9), the result is written 


as 


VI + ^vrrlii »foVV> 

st + p “j iT + 15 “i iT + f^- 1 - 


+ (p Uj ’ U k ’ + p’ Uj « Ufc + p,- ■ u , , } 3ui 


8xk 
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+ (p. Uj/Ufc 1 + p 1 uj/ u k + p f Ui 1 Uk f ) JiL 


+ [(pu k ' + p’u k + p r u k * ? ) a( “l 1u J l) l p „ a( PiV> 

3 x k J 3 x k 

► 2 L’ u. T 3 < P u k + P f u k + P' V k >1 

L 3 3 XT, J 


li-i * + ui* ^ ¥1 

3x£ 9xj _ 


+ ui 


TS rlk 1 

3x k. 


P u k 


: J + 2 u i u i 5^— 


, . 3(pu k ’) , 3(P ui'u^uk') 

u • u-j' — — + J — 

1 J 3 x k 3 x k 

f Diffusion of the correlation u^’uj' 

by the turbulent motion itself 

This term is dropped for the 
incompressible flow, since 


3 (uk + u k ’) 


0 , (A) 


9(uk + uk*) _ q (B) 

3 x k 


(A) ' (Br ^ = ° 
dX k 


- 83 - 



Substitute E^a (B-6) into j we obtein 


u k 


acvj^J 

9s k J 


u i u j 


9(pu k 

3x k 


3 



|3(p’u^ •}• p'ufe’ - p'uk') 
U H 3xk 



9 (pui'uj ’uk’) 
* 9xk 


ihus iiqo (iJ-iO) is rewritten as follows i 


- ^Cttj’uj * ) 3ui 3uj 3(p’ui ? ui ") 

p at +p ’ u d , at + p,u i* aT + ' §1 


+ 


(pUk) * Uj ' 



+ 


(pUk) f Ui 




(pui'uj ’uk’) 

9 x k 


+ (p'u k + p 1 « k ') 'A 


8xk 


1 


-J- p u k 


aCui’ui’) 

~*r~ 


+ 



a(p*uk + p^uk 9 - p’uk*) 1 f , 

a»k J ~ u± aT 


»(p ? uj > 

3xi 


+ 


a X j 



hll 

3xi 


+ 


t 


3ui' 

axj > 



The general pressure The tendency- 
diffusion terms towards isotropy 

term 


3Tik' 

3xk + U i 


(B-12) 


3xk 


(B-13) 
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Since 


(P 


_ aCuj/uj') 

*u k + P'uk’) a Xk 


*' 9 (p’up 3 (p’u k ’) 

f it j ^ ■ X ii j T it j 


+ Ui ! 


U J 3*,, + U1 U J 3*k. 


3 (p' UkUj'uj*) + SCp'uk’uj’u-j ') 

3 x k 3 xk 


( 8 - 14 ) 


Eq. (B- 13 ) is written for steady flow as 


, .r, r 3ui . ■ r^Ui , d<pul’ui ’Uk’) 

(pu ^ ^ + (pUk) Ui ^k 3 xk ' 


3(p’ukui’uj ') SCp'uf’u^ r Uk') __ 3(ui*u-j') — 3p'uk' 

+ P uk 3 xk - U 1 ,U J’ 9 x k 


+ 3xk + 3x k 



3(p*uj') 

acp'ui') | 

/ 3uj ' 3uj\ 


3rik T 3xjk* " 

T ~ 

3x± 

+ ^ J + 

^ ' 9*i 9xj / 

+ 

|_“3 9xk + Ui 3 *k 


(B- 15 ) 


where 


acpui'uj 'u k ') a(p , u k u i , uj’) acp'ui’uj 'ujj') 

9 xk + 3 x k + 3 x 1 c 


9 I(pu k ) 'ui’uj '] 
3 x k 

CB- 16 ) 


and 


-ui’uj ' 


3 p’u k ' 3 [ (P 1 ^ 1 ) (Ui’uj ’)] Sui’uj' 


-+ p’u k ’ 


3 xk 


3 xk 


3 x k 


(B- 17 ) 
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Thus Eq. (B-15) is 


(Pu k )'u. 


3u. , 

5.. f x + 




3u, + 3 t (puk) ’ui'uj ' ] 3[ (p’uk’) (ui’uj’)] 




+ (pu k + P’u k ’) J. 




SCp'uj*) 3(p t u i ’) 
3x i ‘ ~ 3s, 


+ P 


V + 3ui ! 



H 3 *1 


3tik' 3rjk’ 


uj “5 + u ’ -r 


(B-18) 


Let us assume 


3[(pu k )»u i »u j ’] 3[(p’u k ? ){u i f u j T )] 


3 si 


: k 


3xi 


k 


3(p Ui'uj ’u k c ) 


(B-19) 


then we have a system of turbulent fluctuation equation in a simple 


form as 



(pu^'u, 8 

+ (pu k ) 'ui' . 


3xk 


3x k 

= _ 

3( F 5 u i f ) "j 

+ ‘ ,/!±' 

L 3 Xi 

dXj J 

* L * 1 ' 


3s k 


u k +PX*) 

3x k 


-)M, 


3Tik’ 3xjk 

"ST***’ 


I ~ 


(B-20) 
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As Eq. (B-2) showed 


Tik = Tik + Tik» 


. u^Hl + ^Hk _ 1 6 ML)+ r 6.. SSL 
y Va*k 3*i 3 °ik 8xjt/ C ik 3x£, 

thus neglecting the bulk viscosity term, we have 

. /3^ J - 1 3uk * 2 9u2' \ 

Tik ’ = + " 3- 6 * k fcj ) 

Let us assume 


(B-21) 


uj* 


3x k 


/ 3u% 
U ( 


_ I fi±k 3u£ 

3 3x2 


1 \h , 3 T ( 3u* 

) ui 3x k V {^ 

L 


3 Xj 

2 3u2 ,; 3_ „ 

’ ‘ 3 ^ Mj ° 


(B-22) 


then the last term oh the left hand side of Eq. (B-20) becomes 


V 


9t 


ik' 


3x 


k 


+ u i’- 


3x 


ik 


3xi 


= u ' 3 , 

L 9u i’ \ 

l+Ui'- 

3 / 
( V 

3x k V 

3 uj 

J 3x k 

r 3x k J 

3x k 

= -3_/ u 

3ui ' uj 1 

)- 2 y 

3ui f 

9u j 

3x \ 

3x k ' 

X 

CO 


3 / 

Sui'uj 1 ’ 

)- 2y 

3ui' 

3u j ' 

3x k V 

3x k y 

3xk 



(B-23) 
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Therefore Eq„ (B-2Q) is simplified to 
[Turbulent Fluctuation Equations] 


(puk) *uj ’ 



(P^k) 


chl-j 

9sk 


9 (p Ui' 


J j u k 


3x k 


(puj+p'uk 1 ) 


&U£ f U-j 

9xk 


l*i P’uj* 9p r ui f \ 

“Syl 


\ 3x i 3*j) 9^k l y 


3*k / 


2 y 


9ui T auj* 

9?k 3xk 


(B-24) 

Add three equations above for i-j- 1, 2s and 3 considering the 

boundary layer assumption, then we obtain the following turbulent 
kinetic energy equation; 


[Turbulent Kinetic Energy Equation] 


-- 9K 


3K 


Pu £ +(pv + p'v') | = ~ 2(pv) 'u' 


9y 


3y 


9K 

3y |/ Sy 


- (pK 0 -J-2P s )v’ 1 - 2 9v'9v^ ^ 3w ; 3w 


\3y 3y ’ 9y 3y 3y 


dw \ 

dyj 


(B-25) 


where 

and 


Kq - -f- V P V ? i 


K ^ u 9 u ^ + v 7 v p ^ w 1 w ? 
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